 &MODEL
 NLIST   =           1,
 NSTEP   =        2000,
 EPSI    =  1.0000000E-07,
 WGHT    =   1.000000    ,
 NOBS    =         120,
 NVAXZ   =           4,
 NVAX    =           4,
 NORDER  =           1,
 NVAZ    =           0,
 NLAMZ   =           0,
 NLAMY   =           1,
 NLAMX   =           1,
 NSCALE  =           1,
 NEXPY   =           0,
 NRES    =           0,
 NCOR    =           0,
 NCOV    =           0,
 VARNAME = LOG+AU1  ,
 NVARIANT        =           2
 /
 &AUTOCOR
 RHO     = 2*0.000000000000000E+000  ,
 PI      =   1.00000000000000     ,  0.000000000000000E+000,
 NTPI    =           0,   538976288,
 MATFILE = R320km.in                                                                                                                 
                                       
 /
 &YBOXCOX
 VLY     =  0.000000000000000E+000,
 NTLY    =           0,
 NUMLXY  =           0
 /
 &XBOXCOX
 NUMX    =           1,           2,           3,           4, 46*0,
 NUMLX   = 4*1, 46*0,
 VLX     = 10*0.000000000000000E+000  ,
 NTLX    = 10*0
 /
 &YXZDATA
 YXZFILE = Canada1976.in                                                                   ,
 YXZFMT  = (5e15.5)                                                                        
 /

VARIABLE       MINIMUM         MEAN          MAXIMUM   NULL OBS.      SCALE
DEPENDENT
TOTFLOW     0.2577000D+04  0.1476588D+06  0.2736000D+07     0     0.1000000D-05
INDEPENDENT
POP         0.8591300D+05  0.5297039D+06  0.2802800D+07     0     0.1000000D-05
REVE        0.3698300D+04  0.4258776D+04  0.4808200D+04     0     0.1000000D-02
LANGFR      0.5600000D+01  0.7329333D+02  0.9990000D+02     0     0.1000000D+00
UTILITY     0.5655600D-04  0.4915871D-01  0.1172800D+01     0     0.1000000D+01
CONSTANT    0.1000000D+01  0.1000000D+01  0.1000000D+01     0     0.1000000D+01


MRS of X(l) with respect to X(k) - at the sample means
===================================================================================================================================
        X(k) POP         REVE        LANGFR      UTILITY     CONSTANT  
X(l)
===================================================================================================================================

POP          0.1000D+01  0.7451D+02  0.3456D+03  0.4914D+07  0.1889D+06

REVE         0.1342D-01  0.1000D+01  0.4639D+01  0.6595D+05  0.2536D+04

LANGFR       0.2893D-02  0.2156D+00  0.1000D+01  0.1422D+05  0.5466D+03

UTILITY      0.2035D-06  0.1516D-04  0.7034D-04  0.1000D+01  0.3845D-01

CONSTANT     0.5293D-05  0.3944D-03  0.1829D-02  0.2601D+02  0.1000D+01


INITIAL ESTIMATES

             LAM    PARAMETER      STD-ERROR        STUDENT-T    GRADIENT       DERIV.Y         ELAST.Y
                                                                                DERIV.E(Y)      ELAST.E(Y)


POP           X1  0.1389205D+01  0.5457770D-01     25.4537107  0.4002112D-09  0.3872510D+00      1.3892049
                                                                              0.4970914D+00      1.3892049

REVE          X1  0.8322593D+00  0.6932686D+00      1.2004861 -0.6521993D-09  0.2885581D+02      0.8322593
                                                                              0.3626500D-01      0.8322593

LANGFR        X1  0.6643968D-01  0.5550501D-01      1.1970033 -0.8357800D-09  0.1338513D+03      0.0664397
                                                                              0.1916400D-02      0.0664397

UTILITY       X1  0.6335200D+00  0.1496708D-01     42.3275532  0.2744231D-08  0.1902915D+07      0.6335200
                                                                              0.2843971D+02      0.6335200

CONSTANT      X0  0.4955077D+00  0.9808520D+00      0.5051808 -0.4526935D-09  0.7316608D+05      0.4955077
                                                                              0.9221363D-01      0.4955077

VAR(W)            0.1312324D+00  0.1694202D-01      7.7459667  0.0000000D+00

INITIAL VALUES
 
RHO1           0.0000000D+00
PI1            0.1000000D+01
 
LAM(Y)         0.0000000D+00
LAM(X1)        0.0000000D+00


LOG-LIKELIHOOD =  0.3396262D+03

NUMBER OF OBSERVATIONS              =            120

PEARSON-R2(E) (unadjusted)          =         0.8850
PEARSON-R2(E) (adjusted for D.F.)   =         0.8810

PSEUDO-R2(L)  (unadjusted)          =         0.9441
PSEUDO-R2(L)  (adjusted for D.F.)   =         0.9422

PROBABILITY of Y to be at the limit =         0.0000
MEAN VALUE of E(Y)                  =  0.1860993D+00

MEAN of residual W                  = -0.4950736D-12
VARIANCE of residual W              =  0.1312324D+00
95% CONFIDENCE INTERVAL for residual W
Lower Bound = -0.7100297D+00   Upper Bound =  0.7100297D+00




initialization cards for control routine follow - -  id card - no. parameters , no. steps , step size 
minimization tolerance - starting values of parameters - relative uncertainties of parameters (error estimates) -  
and starting direction  (relative parameter increments for 1st step)
identification line
         1      2000   0.1E+01   0.1E-06
  0.00000E+00
  0.10000E+01
  0.00000E+00

first entry to fcn , parameters of current step follow - - f,x(i)
  -0.3396262428D+03  0.00000D+00

new absolute minimum , parameters for fcn entry     6  follow - - f,x(i)
  -0.3632829946D+03  0.96598D+00

new absolute minimum , parameters for fcn entry     7  follow - - f,x(i)
  -0.3632829946D+03  0.96598D+00

step to unallowed region , parameters for fcn entry    12 follow - - x(i)
                     0.15986D+08

step to unallowed region , parameters for fcn entry    13 follow - - x(i)
                     0.79932D+07

step to unallowed region , parameters for fcn entry    14 follow - - x(i)
                     0.39966D+07

step to unallowed region , parameters for fcn entry    15 follow - - x(i)
                     0.19983D+07


the run has been terminated since no improvement can be found on the current minimum,
 the parameters of which follow   - - f,x(i)
  -0.3632829946D+03  0.96598D+00


all numbers of last entry refer to final pt. of run


the lowest value of the function (fcn) was  -0.3632829946D+03 for entry   24 - - x(i) follows
  0.96598E+00



run statistics follow - -

 nfcn=   25 nst=    5


MRS of X(l) with respect to X(k) - at the sample means
===================================================================================================================================
        X(k) POP         REVE        LANGFR      UTILITY     CONSTANT  
X(l)
===================================================================================================================================

POP          0.1000D+01  0.1577D+03  0.1700D+04  0.4335D+07 -0.7709D+07

REVE         0.6340D-02  0.1000D+01  0.1077D+02  0.2748D+05 -0.4887D+05

LANGFR       0.5884D-03  0.9281D-01  0.1000D+01  0.2551D+04 -0.4536D+04

UTILITY      0.2307D-06  0.3639D-04  0.3920D-03  0.1000D+01 -0.1778D+01

CONSTANT    -0.1297D-06 -0.2046D-04 -0.2204D-03 -0.5623D+00  0.1000D+01


FINAL ESTIMATES

             LAM    PARAMETER      STD-ERROR        STUDENT-T      STUDENT-T    GRADIENT       DERIV.Y         ELAST.Y
                                   (CONDIT.)        (CONDIT.)      (CONDIT.)                   DERIV.E(Y)      ELAST.E(Y)
                                                    (PAR = 0)      (PAR = 1)

POP           X1  0.1418847D+01  0.5059234D-01     28.0446973                 0.2594500D-08  0.3955139D+00      1.4188468
                                 0.5059234D-01     28.0446973                                0.4505784D+00      1.4188468

REVE          X1  0.1799387D+01  0.7544361D+00      2.3850752                 0.1651441D-08  0.6238772D+02      1.7993868
                                 0.7544361D+00      2.3850752                                0.7269231D+02      1.7993868

LANGFR        X1  0.3336479D+00  0.8172845D-01      4.0823961                 0.8454245D-09  0.6721765D+03      0.3336479
                                 0.8172845D-01      4.0823961                                0.9141665D+03      0.3336479

UTILITY       X1  0.5708005D+00  0.2277218D-01     25.0656951                -0.1362232D-08  0.1714523D+07      0.5708005
                                 0.2277218D-01     25.0656951                                0.2623274D+08      0.5708005

CONSTANT      X0 -0.2065001D+02  0.6171206D+01     -3.3461865                 0.1898969D-09 -0.3049155D+07    -20.6500063
                                 0.6171206D+01     -3.3461865                               -0.3571227D+07    -20.6500063

VAR(W)            0.7274127D-01  0.1093632D-01      6.6513496                 0.0000000D+00
                                 0.1093632D-01      6.6513496
 
RHO1      (E)     0.7469335D+00  0.6788673D-01     11.0026443                 0.1993585D-06
                                 0.6788673D-01     11.0026443
PI1       (F)     0.1000000D+01
 
LAM(Y)    (F)     0.0000000D+00
LAM(X1)   (F)     0.0000000D+00


LOG-LIKELIHOOD = -0.1294578D+04

NUMBER OF OBSERVATIONS              =            120

PEARSON-R2(E) (unadjusted)          =         0.9433
PEARSON-R2(E) (adjusted for D.F.)   =         0.9408

PSEUDO-R2(L)  (unadjusted)          =         0.9797
PSEUDO-R2(L)  (adjusted for D.F.)   =         0.9788

PROBABILITY of Y to be at the limit =         0.0000
MEAN VALUE of E(Y)                  =  0.1729407D+06

MEAN of residual W                  = -0.2064958D-01
VARIANCE of residual W              =  0.7274127D-01
95% CONFIDENCE INTERVAL for residual W
Lower Bound = -0.5492732D+00   Upper Bound =  0.5079740D+00


CPU TIME:      0.39 sec
